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Relativistic Astrophysics 
Final Report 
I. Introduction 

This report follows two independent lines of research. The first is 
an investigation into the coordinate independence of gravitational radiation 
discussed in Section II and the second is an extension of the parameterized 
post Newtonian (PPN) approximation given in Section III. We conclude the 
report by speculating on the direction of future research indicated by the 
work completed here. 

II. Gravitational Radiation in Asymptotic de Sitter Space 
Although an investigation of gravitational radiation seems somewhat 
remote from the parameterized post Newtonian (PPN) analysis of gravitational 
theories, its importance as a definitive test between theories is just now 
being recognized.^ Part of the reason for this budding importance comes 
from a partial failure at the level of the PPN approximation to distinguish 

9 

between alternate and viable theories of gravitation. Attempts to firmly 

root Einstein's general theory to experimental results has resulted in a 

plethora of theories that agree experimentally in the PPN approximation. 

Deviations between some of the theories are only expected in the "post" PPN 
3 

approximation. It is at this level of approximation that one begins to 
detect varying contributions to gravitational radiation.^ 
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In general the PPN approximation is an expansion of the metric 
(Lagrangian) to fourth order in the velocity which is small in the neighbor- 
hood of the solar system. In fact near the sun one expects that the 
quantities 

p(x) ~ v 2 ~ U(x) ~ 

where p is the density, U is the potential and p is the pressure, are all 
of the same order of magnitude. Far from the solar system, we expect the 
metric to take its Minkowskiian form: g = (-1,1,1,!). [Greek symbols 

r v 

represent the space-time components 0,1, 2, 3 whereas latin symbols represent 
spatial components 1,2,3.] Thus if one carried out a consistent approxi- 
mation scheme to successively higher orders, gravitational radiation effects 
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would occur in the seventh order. Therefore, a clear understanding of 
gravitational radiation seems necessary in order to interpret such an 
expansion. 

A general unsolved problem exists with the understanding of gravitational 
radiation in the interpretation of coordinates in relativistic theories. 
Although the covariance of the equations seem to indicate an independence 
of the choice of coordinates in calculation, it does make a difference when 
one obtains the "physical components" necessary for comparison with experi- 

r 

rnent. 3 For instance, the formulation of gravitational radiation in one 
frame of reference (theory) may parrot the properties of gravitational 
radiation but in the end, fail the test of coordinate independence. This 
might occur if, for instance, a spacelike coordinate could become timelike 
and vice versa as inside the horizon of a blackhole.^ 
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Here we have investigated the gravitational radiation conditions in 
an asymptotic de Sitter space, a theory conformally equivalent to general 
relativity. Physically, this is an interesting problem for at least two 
reasons: 


(1) The results will help demonstrate the coordinate independence 
of gravitational radiation, i.e., a physical interpretation 
useful in the detection of gravitational radiation. 

(2) It leads directly to an investigation of the group structure 
imposed on space-time by an asymptotic de Sitter universe. 

The second case has its importance in the possible classification schemes 
for elementary particles. ^ The former is a necessary step towards correctly 

o 

quantizing the gravitational field. 

Gravitational radiation in asymptotic spaces has been investigated by 

q 

many, notably the works of Bondi, Van der Burg and Metzner (BVM) in an 
empty, asymptotically flat space (i.e. asymptotically Minkowskiian) and of 
Hawking in a dust-filled, asymptotic conformally flat Friedmann universe 
with negative curvature.^® Here we investigate a matter-less but asymptotic 
de Sitter universe. For completeness, we include the effect of the cosmo- 
logical constant. We use the method of the BVM empty space approach which 
is well suited for this type of investigation. 

To use the BVM method, it is necessary to transform their metric into 
a de Sitter-space metric. This can be done thru a conformal transformation 
of the metric^ ^ 




(1) 
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where 


a = £n 


4R 


( 2 ) 


4R 2 -u 2 -2ur 

where R is the “radius" of de Sitter space, u = t - r is the null coordinate 

and r is the ordinary radial coordinate. Since the BVM metric is axially 

symmetric, so will the de Sitter space metric be axially symmetric. After 
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the correction of several misprints in the work of Bondi, et al^, one 
obtains the main field equations 
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For completeness we list the consistency field equation 


V-n ^11^ ^1^1 ^1^ 
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The origin of the names for the various field equations comes from the way 
in which a solution for the functions y>$>U and V occur. An expansion is 
assumed for y at some instant in retarded time u. Eqs. (3) - (5) then yield 
solutions for 3, U and V at the same time. Eq. (6) then gives the time 
development of y from which 3, U and V can be found for all time. Eqs. (7) 
and (8) then represent relations between the expansion parameters and 
initial conditions (integration constants), and finally Eq. (11) must be 
trivially satisfied providing the solution for y» 3» U and V is correct. 

Since we are working in a cosmological space, in particular de Sitter 

13 

space, the field equations take the generalized form 

R’-lgiR + ^¥ = f. (12) 

yv 2 y yv y yv y\> 

where \ is the cosmological constant. 


In empty space, they become 




( 13 ) 


so that, in general, the empty space equations are not set equal to zero. 
In order to see how the field equations are solved, we note that the metric 
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and then 


Also 


°11 " a l 
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I + ra, — s — <p 


'1 B 


(18) 
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Then (16) now becomes 
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In terms of the variables 
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(20) finally becomes (after dropping primes) 
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where derivatives are with respect to q except for the function 


„ _ 3q _ B-u J. 
q l = 9? ' ~B~ * 


( 23 ) 


which occurs in the charge of variables 
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After using (23) in (22), we get finally 


Q _ 1 „ 2 
e l " 2 q Y 1 
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Since = Xg^ 2 = 0, then in terms of the variables q, u, (4) now becomes 
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Since 


R p 2 (8-y). d e 
r%22 c T «>oo 


^ - - 2q 2 e 2B A 
33 .2 cq e \ 


sin 0 


(26) 
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Since 


then (6) becomes 
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i. 3 i /\„28 


- qW^ = | (^) - 4qa 0 2 B r 4a 0 B 0 


,B-u 2 ,( <+X 2(4' O V)(qB 0 ) 1 Bff’vMq 3 .^) 

+ B >r 2 + ~ 13 + 12 + T 

I Q Q 0 n 


* u ( O )lC0t e ♦ ? (Ct), 


-4|l\ / 4.3 


2 ,* 3 V ( . U(4 j V 2 ), 


(q 2 B 2 U),(q 3 V) (3. Uq> 3 V)p B,Uif 3 V) 

+ — V — + ^ir- lt -^^cote 


3i> \ / . 3. 


2y 1 U(4> j V 2 ) b _ u 2 (c^VK^V)^ 8-j (4> 3 V)(4» 3 V) 1 


- (^)C jr 

B 2q^ 


(qBr)i(^V) 2 


2 y ,,2 



15 


Since R 
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♦ |- ( $ J1 - + (y 1 - b ,)uu 1 
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(33) 


The complex of Eqs. (25), (26) and (28) for g, U, and V respectively, 
can be solved if we assume a power series expansion in 1/q for y of the form 


f(u,e) d(u,e) g(u,e) k(u,e) 


Y = 


+ —4 — + e(q ) 
q 


(34) 


The d term is dropped since it gives rise to a log term in the solution for 
U. The solution for (25) then becomes 


1 f 2 3 - ' 4 

/ x 4 ^ 4 ^9 5 fk 
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q q q 


+ 0(q“ D ) 
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where H is an integration constant. The H term is also dropped since it 
gives rise to a log term in the solution for V. However, this term can be 

g 

reduced to zero by a suitable coordinate transformation. Using (34) and 
(35) in (26), we obtain 


(Y->U ~ L - (f 2 +2f cot 0) ^ + (2N+3ff 2 +4f 2 cot 0) ~ 
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where L and N are integration constants. 


Using (36), the solution for V from (27) is 
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b q Z . 

3 

If we replace the parameter f -> C-c /6, the solution represented by Eqs. 
(34), (35), (36) and (37) will then reduce to the assymptotic flat-space 
results in the limit B -+ °° and x -* 0 providing we put L = o. In asymptotic 
flat-space, the constant of integration L in the equation for U (Eq. (36)) 
must vanish in order to preserve the signature of the metric. In flat-space 
V ~ r and U ~ L. Thus 


Ve 2 ^ 2 2 y 

g = U re q 

y oo r 


would eventually change sign for large enough r. But in asymptotic de Sitter 

q 3 

space, 4) V - 0(q ) and U ~ L; so that g rtrt does not change sign for large q. 
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Thus we cannot a priori set L = 0 as in the flat-space case. The solution 
represented by y, 0, U and V (Eqs. (34) - (37)) are then substituted into 
the time development Eq. (30). This yields the following conditions on 
the expansion: 


and 


4( 


| fX = L 2 - L cot 0 
2 

5jp-)g = 2fM-N 2 +N cot 6 + j kX 


(39) 


- 6g!_2-4g2L-6gL cot 0 


(40) 


With these conditions, the form of y is preserved and the development of 
the system is fully determined from initial conditions provided the functions 
f, k, N, M, L are known. The consistency field equation (33) is trivially 
satisfied to 0(q~^) by this solution for y, 0, U, and V. 

This solution, along with the constraints (39) and (40)., reduces the 
supplementary field equations (31) and (32) to inverse-square form. Setting 
these 0(q ) terms equal to zero yields repsectively the following relations 


- «T-f f o‘ + «T-)Ii < f 22o +3f 2o cot 6 - 2f o> 


-B-u \2j;2./B-u \rl 
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- 2 fg 2 COt 0 +2f2+4f2cot^0+2fcot 0 +2fcot^0) 

- \ L 2 (3M-4f22-10f 2 cot 0 -2fcot 2 0+4f) 

- \ L 2 B2ff 22 +2f 2 2 + | ff 2 cot 0 +f 2 + |Mcqt 2 0> 

- \ LL 2 (8N+ ™ ff 2 + ^ f 2 cot 0)- | L 2 2 f 2 + f L 22 f 2 (41) 
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- X(g 2 +g cot 0)+L 2 (2N+ ff 2 + | f 2 cot 0) 

+ L(3N 2 +3N cot 0 +3ff22+f 2 2+ f ff 2 cot 0 

- 3f 2 - | f 2 cot 2 0) (42) 

where (42) was used in Eq. (31) in order to remove the N q dependence from 
(41). Thus, the time development of M and N are known provided the functions 
g, L, M, N are given for one value of u and f is given as a function of 
u and 6. Thus, we are at the peculiar point that we need to know f, k 
and L as functions Of u and 0. The reason we need to know L is that the 
constraint between f and L from Eq. (39) still leaves L unknown up to an 
arbitrary function of u. We can, however, obtain an equation for k Q 
provided we carry out the expansion of the time development field equation 

O' 

(30) to 0(q ). But this would still leave L arbitrary. 


An alternate approach to this problem is to require that the solution 


has a definite form in the limit of q This would be, of course, the 

asymptotic'de Sitter space form of the metric 


t 

(B-u 2 )' 


B 


,2 _ 4u „2 , 4uB _ , 
g oo * /n .2x2 q .2x2 q , n . 2x2 


(B-u ) (B-u )‘ 


(43 


But from the solution (34) - (37) and the form of g Q0 from the metric (14), 


we have 


B 


9 °° (B-u 2 ) 2 


Ml. X . l 2 3 q 2 


+ + L 2 + L cot 0 - 2fL 2 ]q +[1 + | f 2 X 


- 2f 2 (B+u 2 ) _ 2f 2 L 2 + 2L(f 2 +f cot 0)] + 0(1) J 


Thus in the limit of large q, we must have 


4u 2 _ B 2 r 4(B+u 2 ) X ,2 t 

9 9 L 9 “o“ L J 


(B-u 2 ) 2 (B-u 2 ) 2 B 2 


or 


X 4 = . 2 
3 “ B ~ L 


and similarly 


2fL = L 0 + L cot 0 


(44) 


(45) 
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| f _ 2f 2 L 2 +2L(f 2 + fcot 0)= 0 (48) 

Eq. (46) immediately implies that L is independent of the coordinates and 
is therefore an invariant. But this implies with Eq. (47) that either L=0 
or f a cot 0. The latter case does not have the correct regularity 

9 

properties for y as 0 -+ 0. Thus 


L = 0 



) 


/ 


(49) 


O 

Finally Eq. (48) implies that either f = 0 or B = u /6. Previously we ruled 
out f a cot 6 on the basis of regularity for y. But from Eq. (39) we must 
now conclude that either f = 0 or X = 0. The consistent choice is 


f = 0 


(50) 


Thus by comparison with the asymptotic flat-space case, the character of 
the solution is very different in asymptotic de Sitter space. Indeed, the 
mass aspect (Eq, (41)) now depends on a different "news" function that becomes, 
due to (40), (49) and (50) 


fl = V- (No+ N cot 9) 

0 (B-u 2 ) 2 


24 „ , 1-2 M n , 96 

D 9 a O ^ cot 0 a 

8 0 B-u 2 B(B-u 2 ) 
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Similarly for N Q , we find 


BM 2 4 

N = — s- + o (g 2 + g cot e) 

0 3(bV) (B-u ) 1 


(52) 


9 14 

On the other hand, the radiation condition * in asymptotic flat space 
expressed in (q,u) coordinates becomes 


k(rr') -!*<<*> (^) -^ 


9r 


39 


(53) 


where y' is a function of 1/r; so in the limit of q «>, the right hand 

side will not vanish unless f - 0. Thus the condition for radiation in 

q-coordinates seems to be satisfied if f = 0, and that the description of 

mass or mass-loss now depends consistently on g Q in the supplementary 

conditions, i.e. Eq. (51). This different dependence of y is not surprising 

since the asymptotic space cannot be described by an empty space axially 

symmetric static metric and therefore related to Weyl's form from which 
g 

Bondi, et al . could identify the mass function. 

We conclude that radiation will occur in an asymptotic de Sitter space 
but that the nature of the solution requires that the 1/q dependence of the 
solution for y begins with the 0(q ) term. This implies that in the 
limit B -* °° , we do not obtain the radiative flat-space results. That is, 
the requirement for radiation in asymptotic de Sitter space seems to be more 
restrictive than in asymptotic flat-space. We now are ready to return to 
the argument following Eq. (49) in which we concluded that f > 0. This 
choice was, however, only consistent with a solution in asymptotic de Sitter 
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space. Instead, if we had chosen X = 0, then B •+ °° which would then 
automatically restrict our solution to flat space. Thus, we are left with 
the unusual consequences of radiation in de Sitter space or no de Sitter 
space solution at all. 

III. Gravitational Theories with Non-Zero 
Divergence of the Energy-Momentum Tensor 

Until recently it was generally believed that any viable (acceptable) 

theory of gravitation must be a metric theory. In particular, any theory 

1 C 1 c 

must satisfy the following two metric postulates: ’ 

(1) There exists a metric of signature 2, which governs the proper- 

2 ct B 

length and proper-time measures: ds = g^dx dx. 

(2) The response of stressed matter is embodied in the divergenceless 
energy-momentum tensor. 

1 5 

Based upon this assumption, the Caltech program has vigorously compiled 

and analyzed various theories which can be tested against experiment 

through the parametrized post-Newtonian (PPN) expansion of the metric.^ 

Unfortunately a number of difficulties have arisen in this program. 

Among other problems is that it is now known that there are theories which 

agree with the experimental values of the PPN parameters and hence the 
18 

Einstein theory, but may disagree theoretically with the Einstein theory 
in the as-yet-experimentally-inaccessible next-order expansions of the 

1 9 

metric. Furthermore, there is at least one example of a nonmetric theory 
which approximates a metric theory to first order and is in agreement with 
the Einstein theory to this same limit. Finally, we have shown that if we 
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replace the second condition for a metric theory (although perhaps not the 

15 

more generalized integral form discussed by Thorne and Will in the partic- 

20 

ular case of a Brans-Dicke gravitational theory, we still obtain a con- 

21 

sistent solution in the post-Newtonian expansion of the metric tensor, 

22 

g As a result, we have published a paper where we have modified the 

otp 

Brans-Dicke theory by assuming that the divergence of the energy-momentum 

tensor is proportional to the covariant derivative of the scalar curvature. 

Rastall argues that the form of the divergence of the energy-momentum tensor 

17 23 

is not ruled out experimentally, at least for the Einstein theory. ’ 

However, no ad hoc additions to the usual Brans-Dicke field equations are 

required as in Rastall 's case or as in the steady-state theories of which 

24 

this is a natural possibility. Very general approaches to scalar-tensor 

25 

theories obtainable from Lagrangians have been considered by Bergmann and 
26 

Harrison. In particular, Harrison shows that many theories thought distinct, 

such as the Einstein theory, Brans-Dicke theory, and some steady-state 

theories (among others), are in fact conformally equivalent. This implicit 

embodiment in the field equations of general relativity of various physically 

inequivalent but conformally equivalent versions of the scalar-tensor 

27 

formalism has been recognized by others. Dicke showed that even a unit 

28 

transformation has this effect on the scalar-tensor theory. McCrea 

recognized that formally the steady-state theory of Hoyle logically has 

the same structure as the Einstein equations if one adds a suitable term to 

T^ v . Indeed, Hoyle's original approach^ to the steady-state theory is 

equivalent to the assumption that T^ v . f 0. However, his later identi- 

»u 

29 ■ 

fication of this source with a scalar creation field does not apply here. 
Formally we differ in that no ad hoc additions to the field equation or T yv 


25 


are required as in Hoyle's steady-state theory or as in Ras tail's case. 

i 

The extra degree of freedom of the scalar field determines uniquely the 

field equations for which we need only supply a reasonable assumption for 

the divergence of such as done by Rastall. It is thus noteworthy that 

the Caltech group is now actively considering nonmetric as well as other 

30 

possible theories of gravity. 

In a second paper, we have discovered the formal equivalence of 

massive Brans-Dicke theories of gravitation with Brans-Dicke theories with 

non-zero divergence of the energy-momentum tensor. As an added bonus, we 

find that this equivalence constrains to < 49. 

In order to demonstrate the general theory we carry out in detail the 

4 • 

general calculation of the PPN parameters for two different modifications, 

as yet unpublished, of the generalized scalar-tensor theories referred to 

32 

as the Bergmann-Wagoner theory. In the first case we let 


T pv = o R ;y 
;v 8ir 


(54) 


where a is an unknown coupling constant. As a matter of distinction, we 
note that this modification immediately implies a correction to the 
continuity equation 


If + v • (pv) = 0 


by adding a term proportional to 


cc 



(55) 


(56) 
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to the right hand side of (55). This can, of course, be obviated by a 

o 

priori requiring that c = e(v ). However, Euler's equation 

..A 

p(v • V)v + pVU + Vp = 0 (57) 

will still be modified even if o is a second order parameter. A modification 
of the form 





;y 


(58) 


where <J> is the scalar field does not suffer the above problems since the 

right hand side of (58) is obviously sixth order since <J>. is either 

j y 

second or third order depending on whether it is spatial or time derivatives, 
respectively. The modifications of the form (58) satisfy the usual Newtonian 
equations of perfect fluid hydrodynamics. Furthermore, in this case the 
coupling parameter a is completely free. 

In order to obtain the PPN parameters, we must work in the so called 

4 

standard gauge so we can compare with the work of Will and others. For us 

33 

this is equivalent to choosing the Chandrasekhar-Nutku gauge: 
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From these it follows that 
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For the above we expanded 
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We assume the energy momentum tensor has the form 


T wv = (Po+Po^P) uV + pg*" 


( 66 ) 


p Q = energy density of rest mass 


P-P, 


7T = 


= specific interval energy density 


p - total energy density 


p = pressure; all quantities are measured in the frame comoving with 
the fluid element, and where the four-velocity is given by 
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dj£ 

UV = dT 


(67) 


In order to obtain an expansion for we require an expansion for u to 
4th order and u 1 to 3rd order: 


u° = 1 + f- - u + 0(v 4 ) 


u 1 = uV = v 1 + e(v 3 ) 


Then to the required order 


( 68 ) 
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T (2)lj - P 0 vV + p6j j 
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where the superscript in parentheses indicate order (p/r )0(v ) 


We now modify the Bergmann-Wagoner field equations 
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(73) 
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and 
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which arises when we assume 


r M - w n 

M v;y 8tt ;v 


(75) 


where a is a constant of 2nd order in v, and where has the general 


form 


Vv = + B <« , ) 4 V;p* 1P + (76) 


From equation (73) 


dU 1 R = „ 8tt r T V + T y i 

R v “ 2 6 v R 4> LT M v \ v J 


or 


where 


R ■ r Pm + V 


T ; T • T = T ^ 

‘M " ‘M y* 4 . “ <t> y 


(77) 


(78) 


From equation (74) we can solve for , and from equation (76) we can find 
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r - + |E q z + - f(4,)* ;p 4. : 


8 tt 
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Equating coefficients from both sides in (82) and (83) 
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Equations (84) - (88) can be solved for A, B, C, D, F giving 
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The field equations are then 
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Rewriting (90), we get 
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If we expand all terms on right hand side with 


♦ =>S) (92) 

where e = + . . . , % is a constant of order of the 

gravitational constant and once again the 0(v N ) is given in the parentheses; 
and further if we write 
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we find 
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Using Equations (69) - (72) in Equations (94) - (97), we find 
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But (101 and (102) imply that 
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where G is the measurable gravitational constant. To second order,Eq. (91) 
becomes 
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Next we find that 
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By comparison with Will's generalized PPN metric, we find 
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where 
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But the observable G now has a 2nd order correction 
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In order to solve this equation, we use the so called superpotential 
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with the solution 
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Again by comparison with Will, we find 
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By comparison with the Brans-Dicke scalar tensor theories, the only 
difference is in the parameters 


and 


B = e BD + A 


h " 3 2BD ‘ A 
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Experimentally we find that 


- 0.46 < A < 0.64 

As an aside we note that the gauge equations are both satisfied providing 
o is second order. 

As a second example, we assume the same Bergmann-Wagoner field equation 
(73) and (74) but with second form of the divergence of the energy-momentum 
tensor given by Eq. (58). 
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This assumption leads to field equations 




+ 


<*>(<!>) 

8tt4> 





+ 


;y;v 

8tt 


and 


0)+l T _ ^uv 

"yv 3+2 w M 1 6tt 


u) l +2g4> 

3+2to 



_2 . 8tt t 
^ ^ 3+2a)((j>) M 


(D 1 (<t>) + 16TTCJ(j) 

3+2o)(4>) 


4>’ p <j> + 

»p 


24> 2 A'(<j>) : 2<t>X(4>) 

3+2u)(<J>) 


( 122 ) 

023) 


where prime denotes differentiation with respect to the scalar field f. We 
then expand each side of the equations as before 
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From which we obtain 



By (124) and (125) 
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By comparison with generalized PPN metric, we get 
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Similarly from Eq. (126) we find 
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We note that Q allows an added degree of freedom in comparison with experi- 
ment. At present there is no experiment which can detect its value. The 
most probable test will come from lunar laser ranging which depends on the 
combination of parameters 3 + * n this case A drops out leaving only 
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an n dependent term. At present the experiment is still consistent with 

zero oscillation in the radial variable and therefore no strict limitation 
37 

on is possible. 


IV. Conclusion 

Obviously these types of modifications can be extended to more 
divergences of T^ v then those described by Eqs. (54) and (58). Such an 
analysis is the subject of another paper which is forthcoming. In this 
paper, we consider the coupling a = a(<j>) and consider the general con- 
sistency of the field equations with the Bianchi identities (which were 
used to obtain Eq. (81)), gauge conditions and the Newtonian limit of the 
perfect fluid equations of hydrodynamics. We note that the examples 
described above represent viable gravitational theories since their PPN 

parameters fall with experimental limits of those parameters as listed by 
32 37 

Ni. ’ We should also note that the technique of modification indicated 
here and applied to scalar-tensor theories is in general applicable to 

30 on 

other theories such as vector-metric theories or double metric theories. 
The degree of freedom from a second field allows unique field equations. 

A more interesting problem is presented by the "new" PPN parameters 
of Will and Nordtvedt^ 0 in which for the usual Berginan-Wagoner theory 

y = 1 , (3 = 1 + A 

o^i ~ (*2 - a 3 - £-j - ^2 ~ C3 “ ~ 

The parameters now represent preferred frame effects whereas non-zero 



indicate non-conservative theories. However, in the modification given 
by Eq. (58), 

y=l ,8 = l+ A + fi 
a.j = 0 , i = 1 , 2 , 3 
Cj = 0 , j = 1, 3, 4 
C 2 = 23 + 2 gg - 3y - 1 
= 4 n 

But this says that these types of theories are non-conservative gravitational 
theories. As far as we know, these theories which appear to globally 
violate conservation of total momentum are the first such theories that 
are also viable alternatives to general relativity. The exact nature of 
this result is the subject of a forthcoming investigation. 
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